Charged particle motion in a time-dependent flux-driven ring: 
an exactly solvable model 
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We consider a charged particle driven by a time-dependent flux threading a quantum ring. The 
dynamics of the charged particle is investigated using classical treatment, Fourier expansion tech- 
nique, time-evolution method, and Lewis- Riesenfeld approach. We have shown that, by properly 
managing the boundary conditions, a time-dependent wave function can be obtained using a general 
non-Hermitian time-dependent invariant, which is a specific linear combination of initial angular- 
momentum and azimuthal- angle operators. It is shown that the linear invariant eigenfunction can 
be realized as a Gaussian-type wave packet with a peak moving along the classical angular trajec- 
tory, while the distribution of the wave packet is determined by the ratio of the coefficient of the 
initial angle to that of the initial canonical angular momentum. From the topologically nontrivial 
nature as well as the classical trajectory and angular momentum, one can determine the dynamical 
motion of the wave packet. It should be noted that the peak position is no longer an expectation 
value of the angle operator, and hence the Ehrenfest theorem is not directly applicable in such a 
topologically nontrivial system. 
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I. INTRODUCTION 

A charged particle driven by a time-dependent pertur- 
bation in a quantum system is a nontrivial fundamental 
issue [l|, 0, [E IH IS IE 0, @] • One can access the charged 
particle wave function by placing it in a quantum ring 
threaded by a time-dependent magnetic flux. The vector 
potential A(i) associated with the time-dependent flux 
$(t) times the charge q leading to a phase shift propor- 
tional to the number of flux quanta penetrating the ring , 
this is known as Aharonov-Bohm (AB) effect [1, EE 
In adiabatic cyclic evolution, Berry [l2[ first discovered 
that there exists a geometric phase. Later on, Aharonov 
and Anandan (AA) removed the adiabatic restriction to 
explore the geometric phase for any cyclic evolution [l| . 
Time-dependent fields are also used to deal with field 
driven Zener tunne ling , in which nonadiabaticity plays a 
crucial role [3, EE El- 

In mesoscopic systems, a number of manifestations of 
the AB effect have been predicted and verified [13, El, El, 
I20I l2ll . l22j . On the other hand, Stern demonstrated that 
the Berry phase affects the particle motion in the ring 
similar to the AB effect, and a time-dependent Berry 
phase induces a motive force (24|. It was found experi- 
mentally [23j that a quantum ring threaded by a static 
magnetic field displays persistent currents oscillating in 
period of <&o = h/q, the ratio of Planck constant and 
charge of a particle. 

In the present work, we consider a noninteracting spin- 
less charged particle moving cyclically in a quantum ring 
in the presence of a time-dependent vector potential. 
Such a particle motion can be described by the time- 
dependent Schrodinger equation 



time-dependent vector potential A(t), given by 
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Here P — egPg is the canonical momentum operator with 
eg being the unit vector along the azimuthal angle 9; 
L = L z = (f x P) 2 is the canonical angular momentum 
operator in the z direction; / = mR 2 is the moment of 
inertia of the particle; A(t) = A(t)eg is the vector poten- 
tial; and R is the radius of the circular ring. This time- 
dependent dynamical problem can be solved by taking 
into account the Fourier expansion, time evolution oper- 
ator, and Lewis- Riesenfeld (LR) method [25l. |26||. 



II. A CLASSICAL TREATMENT 

We first analyze the time-dependent problem in a clas- 
sical manner. The time- varying magnetic flux induces 
an electric field E = Eeg such that E = —dA/dt. The 
charged particle thus obtain a kinematic momentum in- 
crement during the time interval from to t, namely 

A Pc = A(mv) = m[v(t) - v(0)] = -q [A(t) - A(0)] , (3) 

where p c = mv is the kinematic momentum. It should 
be noted that both p c and qA are not conservative quan- 
tities, while from Eq. © we see that the canonical mo- 
mentum P c is a constant of motion: 

P c (t) = mv{t) + qA(t) = mv(0) + qA(0) = P c (0). (4) 

Comparing the two identities in Eq. ©, we see that 
the result of Eq. ^ is equivalent to 



where the Hamiltonian H(t) is induced by an external 



Al c =I[w(t) -w(0)] 
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where l c = (rxp c ) 2 = Iui indicates the kinematic angular 
momentum, w is the angular velocity, $ is the magnetic 
flux threading the ring, and the fact 



$(*) = 2vr RA(t) 



(6) 



has been used. Also, Eq.(0J lead us to obtain the follow- 
ing relations 

L c (t) = l c (t) + ±$(t) = l c (0) + f $(0) = L c (0). (7) 

Z7T Z7T 

These identities imply that the canonical angular mo- 
mentum L c , defined by (r x P c ) z , is also a constant of 
motion. 

Now we define the writhing number as 



$0 ' 
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where $o = h/q is a flux quantum. We also define 

L c = n Q h, l c (t) = n c (t)h, (9) 

then we have 

n c (t) = no - ra*(f). (10) 

All of these n's are reaZ numbers. 

Now the angular position of the driven particle is given 
by 

c {t) =6 Q+ f ^dr 
Jo 1 
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Here 9q indicates the initial azimuthal angle; and u>q = 
w(0) = noh/I stands for the initial angular velocity. Be- 
low we denote the initial kinematic angular momentum 
Iq = l c (0) for simplicity. 

Hereafter we solve the quantum version of the problem, 
i.e., Eq. (1), using three different methods. The classical 
quantities l c (t) and 6 c (t) will also appear in the expres- 
sions of the wave function. Their roles in the quantum 
problem will be further explored. 



III. A FOURIER EXPANSION METHOD 

The simplest method for solving the time-dependent 
flux-driven problem is the Fourier expansion method. 
The first thing about the system we discuss is that the 
wave function satisfies the periodic boundary condition: 



ip(6,t) =<ip(6 + 2<K,i). 



(12) 



The most general form of tp for the present problem is 
thus written as 



(13) 



where the c n 's are appropriate coefficients to be deter- 
mined by the initial and the boundary conditions. 

Substituting Eq. (fT3"|) into Eq. (TT]), we can find the 
identity 
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Solving Eq. (|14p , after some procedures we obtain 
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/„(*) = exp { -— ^ [nh - qRA{t)) dt\ , ( 1.5 ) 



and thus 



W,t)= 2 c„exp|-^ / [ 
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n — n$ (r)j dr + inO 

(16) 

As a simple example, let us choose 

c n = TV exp [-cr 2 (n - n ) 2 - i6 (n - n )] , (17) 

where N indicates an appropriate normalization con- 
stant; cr, #o a- n d Ho are real numbers, and n is an integer. 
Substituting Eq.(fTf|l into Eq.fle]). we get 



tp(9,t) = Nexp 
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dr + ino9 c (t) 
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where T = 21a 2 /%. Applying the Poisson summation 
formula 
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on the function 
f(x) = exp 



it 
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we can obtain an alternative expression 
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We note that Eq. (jT6|) is the general solution of the 
problem, whereas Eq. (fT8|) and Eq.(j2"Tj) are two different 
expressions for a special solution defined by the c„ co- 
efficients of Eq. ([TTD. When a 2 t/T < 1, it should be 
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noted that Eq. (fT8"l) converges slowly while Eq. (|2Tj) con- 
verges quickly. This means that in the short time limit 
t < T/a 2 , the wave function is better described by a 
circulating wave packet. However, for the case of long 
time limit t ^> T/a 2 , Eq. (fl~8|) has fast convergency this 
is because, in this expression, only the n rs no terms are 
important. If we further assume no being an integer, then 
at large t the wave function approaches to a circulating 
plane wave that characterized by no, namely 



ip(9,t) « iVexp 



* J2 



21 



dr + inoO 



From these findings we conclude that Eq. (|2"Tj) describes 
the short time behavior and Eq. lfTSf describes the long 
time behavior of the ring system when the wave function 
is defined by Eq. (fT7|) . 



IV. A TIME EVOLUTION METHOD 

In this section, we shall present how to get the general 
solution shown in the previous section in terms of the 
time evolution operator U(t). The state \if}(t)) is con- 
nected with the initial state |^>(0)) through 



W)) = u(t)MQ))- 



(22) 



and the wave function il>{6, t) is given by 

i>(9,t) = (6\U(t)\i>(0)), (23) 

where \9) is the 0-eigenket in the Schrodinger picture that 
will be explained later. 

To begin with, we introduce the canonical commutator 

[0(O),L(O)]=ift. (24) 

From this identity, we have 

[9(t), L(t)] = UHt)[6(0), L(fl)]U(t) = ih. (25) 

Utilizing Eq. ([25|) we can derive 

dL(t) [L(t),H(t)] 



dt 



if l 



0. 



(26) 



we thus obtain the identity L(t) = L(0). Following similar 
procedure it is easy to obtain 



d9(t) [9(t),H(t)] L(0) - n$(t)h 



dt ih 
which gives us 



(27) 



9(t) = 9(0) + ^-f ^^dr. (28) 
1 Jo * 

Here we see that the canonical angular momentum is a 
constant of motion. This is consistent with the classical 
results discussed in Sec. ITT] 



From the above results we have 

[H(t),H(t')]=0 



(29) 



for any two times t and t' . Hence the time evolution 
operator is simply given by 



U(t) = exp 



- / H(r)dr 



exp 




n$(r) dr 



(30) 



To proceed further, we define \n) as the eigenket of 
L(0) obeying 

L(0)|n) = nh\n). (31) 
Then we assume that \9) is an eigenket of e l6 ^ obeying 



9 (°)im = 



\9)=e*\6). 



(32) 



The orthogonal conditions of the two eigenkets can thus 
be expressed by 



1 
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Am-n)e de = Smn (33) 



= 5(6l-6i'). (34) 



These two orthogonal conditions can be derived from the 
closure relations 
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and taking into account the definition 



(9\n) = 



1 



/27T 



=e = n 



(36) 



It should be noted that both 9 and 6' are defined in 
the interval [0,27r). In the coordinate representation, the 
e is an eigenfunction of L lcp = —iUd/d9 with corre- 
sponding eigenvalue nh. This result can be expressed 
as 



(6\L(O)\n)=L Tgp (0\n) = nh(0\n). 
The wave function tp now can be calculated: 

OO 

= E <f|^(*)|n><n|^(0)> 
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If we define 



2tt 



(39) 



then the result of Eg. ([55)1 becomes that of Eq. (TT6")) . 

Using the time evolution operator U(t), we have in- 
deed found the general solution of Eq. (fT6")) . Based on the 
commutativity of the Hamiltonian operator at different 
times (Eq. ([29|) ). the U(t) operator can be constructed 
straightforwardly by simple integration. 



V. THE LEWIS-RIESENFELD METHOD 

In this section, we briefly review the LR method and 
then apply it to solve the present problem. We shall show 
that LR method is not directly applicable, however, a 
simple modification concerning about the boundary con- 
dition makes it applicable to solving the problems with 
periodic boundary condition. 

Traditionally, to utilize the LR method [26j solving a 
time-dependent system, we have to find an operator Q(t) 
such that 

ih^=ih^ + [Q,H]=0,, (40) 
and then find its eigenfunction (px(9,t) satisfying 

Q(t)<px{9,t) = \<px(9,t), (41) 

with A being the corresponding eigenvalue. A wave func- 
tion ipx(8,t) satisfying Eq.JTJ is then obtained via the 
relation 

4> x (0,t) = e***® <p x (6,t), (42) 
where a(t) is a function of time only, satisfying 

a x = <p x -\id/dt-H/h)< Px . (43) 
A general solution ij; of Eq. (JTJ is then given by 

V(0,t)=£>( A )MM)' (44) 

A 

where <?(A) is a weight function for A. 

To proceed let us assume the time-dependent invariant 
operator Q(t) takes the linear form [5j, [6j 



Q(t) = a(t)L + b(t)9 + c(t). 



(45) 



in which a(t), b(t), and c(t) are time-dependent c-numbcr 
functions to be determined. 

Substituting Eq. ([IS"]) into Eq. (|4*D|) and solving these op- 
erator equations, we get 



a(t) = a Q 



b t 



b(t) = b , 



(46) 



c(t) 



C() 



n<s>(r)h 



dr, 



(47) 



where do, bo, and cq are arbitrary complex con- 
stants. Furthermore, substituting Eqs. (|46|) and ((47} into 
Eq.gS]), we find 



Q(t) = a L(0) + b 8(0) + c = Q(0). 



(48) 



In other words, the invariant Q in the Heisenberg picture 
is precisely the linear combination of the initial canonical 
angular momentum L(0) and the initial azimuthal angle 
9(0) with an arbitrary constant Co- Note that in our 
system the L operator is also an invariant. 

It is interesting to ask how the eigenvalue A evolves 
in time. Multiplying the factor e la ^ on both sides of 
Eq.flUJ), we get 



Q(t)^x(9,t) = X^x(e,t). (49) 

Partially differentiating the both sides of Eq. (|4"9")) with 
respect to time and using Eq. (|40|) . we find 

(50) 



X(t) = A(0), 



thus A is a constant. 

To find a solution of Eq. (Q]) , we have to solve Eq. JIT 
first. By solving Eq. (|4*Tj) , we get 



ip\(8,t) = exp 



- ( 



\v(t)9 2 



where 



A - c(t) 
a(t) ' 



us b ° 
a(t) 



Substituting Eq. ([STJ into Eq. (gSJ), we obtain 

[t] 2 (t) +ihv(r)] 



a x (t) = a x (0) 



2Ih 



dr, 



(51) 



(52) 



(53) 



where 



t](t) = h{t) - n$(r)a. (54) 

In the derivation of Eq. (f53")) . we have used the following 
two identities: 

,,2 



v (/i — n^Ti) 



(55) 



Here we see that in general ax(t) is a complex function. 

Although the form of tp\{0,i) = e ia ^cj>x{9, t) is in- 
deed a solution of Eq.(fT]), however, it does not satisfy the 
periodic boundary condition [see Eq. (fl"2"|) ]. This problem 
can be resolved by defining the total wave function tp(9, t) 
as the summation of all ipx(9 + 2mr, t) terms: 

oo 

1>(6,t)= J2 ^x(9 + 2mr,t) 



H exp 

n=— oo 



iotx(t) + -n(t)(9 + 2mr) 
n 



- — v (t)(6 + 2nnf 



(56) 
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It can also be transformed to the equivalent form below 
using the Poisson summation formula: 



and 



\/ 2mv{t) 

OO 

Yl exp 



V Q = -777 = 



exp 



iax(t)+i^6 2 c {t)+in 6 c (t) 



n— — oo 



i7i(n — no)^ 
2v(t) 



in (6-6 c (t)) 



(57) 



From these derivations it should be noted that when us- 
ing the LR method, the boundary conditions have to be 
carefully managed, otherwise one may get an incorrect 
result. 



VI. A COMPARISON TO VARIOUS 
APPROACHES 

In this section we shall show that Eq. (f56"| and (|57|) can 
be cast into the forms of Eq. lf2~lj) and Eq. (p~8|) . respec- 
tively. To proceed further, let us first borrow the param- 
eters no and 9 c (t) from Sec. II in combination with the 
results obtained in Sec. V, we have the simple identity 

a(t)n h + b(t)9 c (t) + c(t) = a n n + b 9 + c . (58) 

Comparing this result with Eq. ([SO")) , we find that they 
are very similar. For simplicity, we define 

\ = a{t)n h + b(t)9 c (t) +c(t), (59) 

in combination with Eq. (|59p . it is easy to obtain 



H(t) = n h + v{t)9 c {t), 

lit) = l c {t) + v{t)6 c {t). 
Further, using the identity 

and the identity of v in Eq. (|55j) . we have 
r?=ll + I±{ v 6l). 
Substituting Eq.([S5| into Eq.([5_3D, we get 



(60) 



(61) 



(62) 



(63) 



e ja A (0) exp 



f_if 



21 2h 



1 _ iM 
x i 



(64) 



In addition, by defining e ta *^ and vq as 



il 
T 



ih 

2^2' 



(66) 



a 



(65) 



we can see clearly that Eqs. (1551) and ([57)1 become exactly 
the same as Eqs. (|21[) and (IT5|) . Hence, we have verified 
that the Fourier transform method, the time evolution 
method, and the Lewis-Riesenfeld method are equiva- 
lent when we choose the coefficients c„ as Eq.fTTJ. This 
restriction is not necessary to find the two equivalent 
general solutions Eq.(TTg]) and Eq.J3S|) obtained by using 
Fourier expansion and time evolution methods, respec- 
tively. It turns out that the Lewis-Riesenfeld method 
seemed to be more restrictive. 

It is now interesting to discuss the physical meanings 
of l c (t) and 9 c (t) we have obtained. Although they are 
originated from the classical treatment, however, in what 
sense do they play a role of dynamic variables in the 
corresponding classical system should be further clarified. 
We would like to bring attention that in the Schrodinger 
picture within coordinate representation, the l c (t) is the 
expectation value of L — qRA(t) = —ihd/d6 — qRA(t). 
However, 6 c (t) is not the expectation value of 6 = 6 with 
respect to the wave function obtained in Eq. (TT5]) , instead, 
it is merely the peak position of the wave packet [see 
Eq.©]. 

In other words, the conventional Ehrenfest Theorem 
is not directly applicable in this topologically nontrivial 
system. This consequence is due to the fact that we are 
not able to distinguish the phase between the angle 9 and 
9+2nir. Hence the 9 operator is not well-defined, only the 
e l9 is a well-defined operator, as has been demonstrated 
in Sec. IV. These facts cause the A losing its meaning as 
an expectation value of the Q operator. 

We finally point out the relationship between the prob- 
lem we have considered here and that we studied in the 
previous work 0] . In that work we studied the motion of 
a charged particle in a one-dimensional space subject to 
a time-varying linear potential. By doing a gauge trans- 
formation as that mentioned in Ref. [27|, the Hamilto- 
nian in Ref. d can be cast into the form of Eq. ((T|) , hence 
the two problems are equivalent if we ignore the differ- 
ence of their topologies. Therefore, by ignoring the factor 
caused by the gauge transformation (which contains only 
a function of time), the circulating wave packet solution, 
Eq. (f2"Tj) , can be viewed as the wave packet solution in a 
one-dimensional system (see Eqs. (24) and (40) in Ref. Q) 
being folded into a ring. That is why in Eq. (|56|) the total 
wave function ip(9, t) can be written as the sum of all the 
ip\(9 + 2nir, t) terms. This folding nature of wave packet 
leads to interferences between different ip\(9 + 2mr,t) 
terms. As a result, the expectation value of 9 operator is 
different from the peak position 9 C of the wave packet. 
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VII. CONCLUDING REMARKS 

In this article, we have studied the problem of a 
charged particle moving in a ring subject to a time- 
dependent flux threading it. After analyzing the prob- 
lem in a classical manner, various approaches includ- 
ing Fourier expansion method, time-evolution method, 
and Lewis-Riesenfeld method are considered and com- 
pared. In the Lewis-Riesenfeld approach, by appropri- 
ately managing the periodic boundary condition of the 
system, a time-dependent wave function can be obtained 
by using a non-Hcrmitian time-dependent linear invari- 
ant. The eigenfunction of the invariant can be realized as 
a Gaussian-type wave packet with the peak moving along 
the classical angular trajectory while the distribution of 
the wave packet is determined by the ratio of the coeffi- 
cient of the initial angle to that of the initial canonical 
angular momentum. In this circular system, we find that 



although the classical trajectory and angular momentum 
can determine the motion of the wave packet; however, 
the peak position is no longer an expectation value of 
the angle operator, and the Ehrcnfcst theorem can not 
be directly applicable. 

Recently, possible schemes of the experimental setup 
to explore the quantum dynamics of a mesoscopic ring 
threaded by a time-dependent magnetic flux have been 
proposed by either capacitively coupling the ring to a 
electronic reservoir [281 ] or applying two shaped time- 
delayed pulses [29|. The quantum dynamics in a time- 
dependent flux-driven ring should be achievable within 
recent fabrication capability. 
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